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Order statistics of periodic, Gaussian noise with l/f a power spectrum is investigated. Using 
simulations and phenomenological arguments, we find three scaling regimes for the average gap 
dk = {xk — Xk+i) between the k th and (k + l) st largest values of the signal. The result dk ~ fc 1 
known for independent, identically distributed variables remains valid for < a < 1. Nontrivial, 
a-dependent scaling exponents dk ~ fc'" -3 ^ 2 emerge for 1 < a < 5 and, finally, a-independent 
scaling, dk ~ k is obtained for a > 5. The spectra of average ordered values Ek = (xi — Xk) ~ k 13 
is also examined. The exponent /3 is derived from the gap scaling as well as by relating Ek to the 
density of near extreme states. Known results for the density of near extreme states combined with 
scaling suggest that /3(a = 2) = 1/2, /3(4) = 3/2, and /3(oo) = 2 are exact values. We also show that 
parallels can be drawn between Ek and the quantum mechanical spectra of a particle in power-law 
potentials. 

PACS numbers: 05.40.-a, 89.75.Da, 68.35. Ct, 05.45.Tp 



I. INTRODUCTION 

Extreme value statistics (EVS) was first developed in 
mathematics [11-0] ■ Its importance was soon recognized 
and emphasized in engineering [4|, [5| , followed by finance 
and environmental problems [a-KLOf. Although applica- 
tions in physics appeared relatively late, they cover a 
wide range of fields including cosmology [ll|, 0, spin 
glasses [lH , random fragmentation [l4[ , percolation [l5| , 
random matrices [IB], and, most actively studied at 
present, interface fluctuations [l7rj24|. 

The extreme value in a batch of data is important, 
but its study makes use of only a small fraction of the 
available information. Accordingly, there have been var- 
ious attempts to extend studies towards near extreme 
characteristics, such as density of states near extremes 
[IE HH , first-passage and return-time statistics [H, [2t| , 
persistence [3(|, and record statistics [3ll - [33| . A natu- 
ral extension (which will be the concern in this paper) 
is to consider not only the extreme, but the sequence 
xi, X2, ■ Xk, ■■■ of the 1 st , 2 nd , k th ,... largest, i.e. ex- 
tract information from the order statistics of the system. 

Order statistics has been much studied in mathemat- 
ics [IE . All relevant quantities are known for inde- 
pendent, identically distributed (i.i.d.) variables, and a 
significant amount is also known for weakly correlated 
systems (29[. In physics, meanwhile, results related to 
order statistics are scarce: the order statistics of the 
brightest galaxies was recently proposed to replace the 
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use of standard candles [35|. Another cosmological ex- 
ample [3(| concerns an inequality satisfied by i.i.d. vari- 
ables for the ratio of the average gap {x\ — x%) and the 
standard deviation o\ of x\. Violation of this inequal- 
ity by the brightness data of galaxies leads to the notion 
that the brightest galaxies in clusters are special, in the 
sense that their brightness cannot simply be statistically 
attributed to the tail end of the luminosity distribution. 
An interesting example of order statistics in statistical 
physics concerns the positions of the k rightmost points 
of branching random walks [37j . The importance of this 
work is that it provides the exact order statistics of a 
system of strongly correlated particles. 

For applications, it is clear that a better understand- 
ing of near extreme properties in correlated system is 
required. Here, we shall make steps in this direction 
by investigating the order statistics of Gaussian signals 
x(t + T) = x(t) of period T with 1/ f a power spectrum. 
Depending on the value of a, such signals correspond to 
well-defined physical processes [white noise (a = 0), 1/f 
noise (a = 1), random walk (a = 2), random acceleration 
(q = 4), etc.]. As a increases, the signal changes from 
uncorrelated (a = 0) to weakly correlated (0 < a < 1), 
and then strongly correlated for 1 < a < oo. Thus, the 
characterization of correlations is straightforward, and 
the distribution of the maximum, x± = maxtx(t) — x(t), 
measured with respect to the time average, x(t), has been 
discussed in the literature [l^, 0, [H)]. The present work 
is an extention of Ref.pil] to order statistics. 

Before summarizing the results, a clarification is in or- 
der. Namely, the signal x(t) is continuous for a > 1 and 
thus, while the meaning of the maximum xi is clear, the 
definition of x-i (the 2 nd largest) etc. may not be obvi- 
ous. The key here is to recognize that, as discussed in 
Sec. II., the signal is determined through a finite number 
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of N Fourier amplitudes. The N independent values of 
the signal can be obtained e.g. by sampling at equal in- 
tervals oiT/N, yielding a discretized signal in which the 
2 nd , 3 rd largest etc. are well-defined. 

Our main result concerns the average gap dk = (xk — 
Xk+i) between the k th and (fc+l) st largest. We find that 
it scales with k as 




for < a < 1 , 
-3)/2 for 1 < a < 5 , 
for 5 < a < oo 



(1) 



The above results are first obtained from simulations of 
the 1// Q signals and then phcnomenological arguments 
are also used to derive the exponents. As can be seen, 
there are three regimes and the scaling exponents match 
at the borderline values a = 1 and 5. It is suspected, 
however, that the power laws at a = 1 and 5 have loga- 
rithmic corrections which are not resolved by the present 
simulations. 

The gap scaling (JTJ implies that the spectrum of the 
average values of the k th largest behaves as 



£fc = \xi 




for < a < 1 , 
for 1 < a < 5 , 
for 5 < a < oo . 



(2) 



The above spectrum for 1 < a < oo can also be ob- 
tained by relating Sk to the density of near extreme states 
which, for periodic signals, is known to be given by the 
distribution of the maxima of the signal relative to the 
initial value (26[. This distribution has been studied in 
detail in (2(| and, in particular, the a = 2, 4, and oo 
cases were solved exactly. Thus, provided the assumption 
about the scaling form e ~ fc' 3 is valid, the exponents for 
a = 2, 4, and oo cases in ([2]) are exact. 

Interestingly, the same spectrum (j2J) emerges in the 
quasi-classical limit of a quantum mechanical system. 
Namely, the energy spectrum of a particle in a one- 
dimensional potential U(z) ~ \z\ e is given by ~ 
£,20/(2+0)^ Thus, there is a one-to-one correspondence 
between the exponents of the order statistics spectrum 
for 1 < a < oo and of the quantum mechanical energy 
spectrum for < 9 < oo. 

In order to arrive at the above results, we begin with 
a short introduction of l/f a signals (Sec. [TTJ) followed 
by the description of the simplest case of a = (i.i.d. 
variables) in Sec. Mil Numerical simulations together 
with scaling considerations for general a are discussed in 
Sec. IIVI The results are rederived from the point of view 
of density of near extreme states in Sec. [V] and the rela- 
tionship to quasi-classical quantum spectra is presented 
in Sec. ED 



II. PERIODIC, GAUSSIAN l/f a SIGNALS 

Detailed discussions of periodic, Gaussian l/f a sig- 
nals and their correlation and roughness properties can 



be found in Refs.[24|, |39(. Here, we just define the rel- 
evant notation and describe how the order statistics is 
evaluated numerically. 

The configurational weight of a Gaussian 1/ f a signal 
x(t) of periodicity T is given by 



V[x{t)} oc e 
where the effective action 

S[{c n };a] = (2Tr) a T 1 - 



-S[x(t)] 



N/2 

71=1 



(3) 



(4) 



is defined through the Fourier coefficients, c„, of the sig- 
nal 



N/2 
n=— JV/2+1 



2ixint/T 



(5) 



Here iV is a positive, even integer, and the phase of c n is 
drawn randomly and uniformly from the interval [0, 271"]. 

From Eqs. © and (|3|) one sees that the amplitudes of 
the Fourier modes are independent, Gaussian distributed 
variables, but only for a = are they identically dis- 
tributed. This is also apparent from the mean square 
amplitude (|c„| 2 ) oc \/n a which is consistent with a l/f a 
power spectrum and independent of n only for a = 0. Al- 
though the Fourier components c„ are uncorrelated, the 
corresponding time signal x(t) is correlated at different 
times t and if for a > 0, and the correlation increases 
with increasing a. In particular, the correlation function 
(x(t')x(t' + t)) is bounded for < a < 1 (weakly cor- 
related regime), while it diverges for a > 1 in the limit 
T — > oo with t/T finite (regime of strong correlations) 

The action in Eqs. © and Q may be formally written 
in the continuum limit as 



S[x(t)} = 





d a / 2 x 


[ dt 


dt a / 2 


'o 



implying the stochastic equation of motion 

d a ' 2 x 



dt a / 2 



m, (atw)) = s(t-t') 



(G) 



(7) 



where is Gaussian white noise with zero mean. In 
this form it is transparent that the signal for a = 0, 2, 
and 4, describes white noise, random walk, and random 
acceleration, respectively. 

It is important to note that since the maximum fre- 
quency appearing in the sum ([5]) is N/T, the series 
docs not resolve fine structure on time scales less than 
t = T/N. Thus, we may view the signal x(t) as a batch 
of N variables x n = x(nT/N) (n = 0, 1, N - 1), and 
determine the order statistics by ordering the x n . As 
usual in extreme statistics, we should consider the limit 
of the batch size going to infinity (TV — > oo) and search 
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FIG. 1: (Color online) Probability distribution of the k = 1, 2, 5, 10 maxima (right to left) for (a) a — 0, (b) a = 0.5, (c) a = 1, 
and (d) a = 2, centered and scaled according to the mean and standard deviation of the k = 1 maximum. TV = 16 384 in all 
cases and the number of simulations (10 6 ) is large enough so that the error bars are of the order of the width of the lines drawn. 
The exact limit distributions are also plotted for x\. One can see that the finite size corrections are large for a < 1. 



for finite results after appropriate rcscalings. In our case, 
the limit of the batch size going to infinity (TV — > oo) is 
equivalent to taking the limit T — > oo with r = T/N kept 
fixed. 



We now describe how the order statistics is evaluated 
numerically. For a given a and TV, the amplitude and 
phase of the Fourier coefficients in ([5]) are sampled from 
their (Gaussian and homogeneous, respectively) distri- 
butions, and the signal x(t) is generated by fast Fourier 
transform. The values x n = x(nT/N) are ordered and 
the largest x%, 2 nd largest X2, etc. are determined. The 
process is repeated 10 6 times in order to obtain the distri- 
bution P/v(a;fe), and to have well defined averages {xk)N, 
as well as gaps dk,N = {xk) n — {xu+ijN- Signals are gen- 
erated in such a way for TV = 16,32, 16384, and the 
TV-dependence of the above quantities is examined. Fi- 
nally, appropriately scaled quantities are introduced, so 
that a finite structure emerges in the N — > oo limit. 



III. ORDER STATISTICS FOR a = 

We begin with the a = limit, with N i.i.d. variables 
drawn from a Gaussian parent distribution. The results 
for order statistics are known [29j in this case. In par- 
ticular, the k th maximum (xk) and its root-mean-squarc 
fluctuation = J — (xk) 2 scale for large TV as 

(x k )~Vh^N , ak~l/VteN. (8) 

The above scalings suggest that, in the TV —> oo limit, 
finite distribution functions are obtained by introduc- 
ing the scaled variable z through x = ajsrz + fow, where 
a at ~ 1/VlnTV and 6 at ~ y/hiN. Indeed, using the "ex- 
perimental" scaling z = (x— (xi))/<Ji, the limit distribu- 
tion Pk(z) for the k th maximum is given by |2!| 

Pk{z) = — — — -tt exp (~kz — e~ z ) , z = az + 7 (9) 
(fc-1)! 

where a = tt/\/6 and 7 is Euler's constant. 
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The Pk(z) distribution functions obtained from simu- 
lations of TV = 16 384 independent modes are displayed 
in FigOJa) for k = 1 [Fisher-Tippett-Gumbel (FTG) 
distribution], k = 2, 5, and 10. As one can see, even 
N = 16 384 is not large enough to arrive at the limit dis- 
tributions for a < 1. This is related to the notoriously 
slow (logarithmic) convergence, as discussed in detail in 
[1113- 
Using @ to express (zk) via (zk+i), the gap between 
the k th and (k + l) st largest is given by 

dk = (zk) - (zk+i) = —r • (10) 
ak 

The spectrum of the average values Zk can be evaluated 
as the sum of the gaps, and one finds for large k 

e k = ( Zl ) - (z k ) =Y. d t = lY.- e ~l lnk - ( n ) 

One should note here that both the distribution Pk(z) 
and the results for the gap and the spectrum remain 
valid for all parent distributions whose tail extends to 
infinity and decay faster than any power law (domain of 
attraction for the FTG universality class). 

In the rest of the paper, we shall consider the case 
a 7^ and concentrate on Pk(z), dk, and e k - In ex- 
plaining the simulation results, we shall make use of a 
phenomenological argument for calculating dk, which we 
demonstrate now for the i.i.d. case. 

Let us assume that we are studying the order statis- 
tics in the unrescaled variable (x\, X2, •■•), and we discover 
from simulations that (xi) ~ vlrjJV and (x& — Xk+i) ~ 
dk/VhiN. Summing over N/2 gaps gives the distance 
from the largest to the to the median of the parent dis- 
tribution [(xjv/2) ~ O(l)], such that 

N/2 N/2 

Qik -aifc+i) w ^dik/VklVw (si) « VkW. (12) 
fc=i fc=i 



Multiplying by y/ln N, the above equalities yield 

N/2 

J^dk^lnN. (13) 
fc=i 

Assuming that the above sum is dominated by the dk 
with simple power- law behavior dk ~ k~ s , one can read 
off the scaling of the gap 

d k ~ k' 1 (14) 

in agreement with the exact expression (jlOl) . 

The scaling assumption, backed up by simulations, can 
thus be used to develop relationships between observed 
quantities. 



IV. ORDER STATISTICS FOR a # 

For a / 0, we shall present the numerical results in 
the following order: First, the probability distribution P k 
(FigO}, second, the TV-dependent scaling of (xi), <J\, and 
(xi — x 2 ) (Fig[2j) are discussed. Then, the average gap 
dk = ( x k — x k+i) as a function of k (Figj3|) is described 
with the iV-scaling removed when necessary. Finally, the 
spectrum e k is obtained by summing up the dk- 

The numerical results indicate that, depending on a, 
there are three distinct regimes of order statistics as we 
shall discuss below. 



A. Weakly correlated regime (0 < a < 1) 

The distributions Pk(z) obtained for a = 0.5 [FigQJb)] 
are characteristic of the results. As can be seen, the 
distributions P k {z) in (b) are rather close to those in (a). 
In fact, the limit distribution of the maximum xi, in the 
interval < a < 1, is known to be FTG distributed just 
as in the i.i.d. case [4l|. At the same time, however, 
the finite-size corrections are large in this regime f24j ]. 
and it is not surprising that the deviations from the limit 
distribution are significant even for N — 16 384. 

FigsQTa),(b) reveal that the widths a k of the distri- 
butions Pfe(x) are of the same order as those for Pi(x) 
and, since it is known that u\ ~ 1/ vlnTV, we have Ofc ~ 
l/Vhi N. Furthermore, the gaps (xfc — Xk+i) are also of 
the same order as o\ and so (x& — x^+i) ~ l/\TnJV. 
This is more precisely checked in Fig|2l where the ratio 
(xi — X2) /(Ti is plotted and, indeed, is of the order 0(1). 
Similar results are obtained for all the (x& — Xk+\}/cr\ 
ratios we examined (k < 10). 




FIG. 2: (Color online) Difference between the k — 1 and k = 2 
average maxima as a function of N for increasing a (top to 
bottom). The asymptotes are indicated on the right. 
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Once it is established that (xk — Xk+i) ~ 1/vlnJV, we 
can plot dk = (xk — Xk+i)V^iN as a function of k to 
examine the scaling properties of the gap dk- As Fig[3] 
shows, dk ~ in agreement with the expectation that 
the weak correlations in the < a < 1 regime do not 
affect the extreme statistics properties of the signal as 
compared to the i.i.d. case. 

It should be noted that the (cc fe — x k+ i) w d fc /VlniV 
scaling in conjunction with the mathematical result 
(xi) ~ VhT/V allows one to repeat the phenomenological 
argument given in Scc lIIII [sec eqs. (fT2|) - (fT4| ] to establish 
that dk ~ A: -1 and ~ In ft (fTTj) . valid in the large ft 
limit. 
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FIG. 3: (Color online) Average unsealed gap dfc for a = 
0, 0.5, 1, . . . , 6 (bottom to top) obtained from simulations with 
N = 16384. The values of the limiting slopes are indicated at 
the end of the lines. The straight line fits are excellent except 
for a = 1 and 5 where logaritmic corrections may be present. 



B. 



Regime of non-trivial scaling (1 < a < 5) 



For a > 1, we enter the strongly correlated regime 
where fluctuations diverge with system size. The large- 
fluctuation regime can be further divided according to 
whether the second derivative of the signal (playing an 
important role in determining the order statistics) is con- 
tinuous (5 < a < oo ) or not (1 < a < 5). Sec II V Bl will 
be devoted to 1 < a < 5. 

We begin with the borderline case (a = 1) separating 
the weakly and strongly correlated regimes. The his- 
tograms are plotted in FigQJc) and, compared to a < 1, 
they draw closer to each other both in terms of location 
and scale. Since the a = 1 case lies at the threshold 
between the two different scaling regimes for a < 1 and 



a > 1, the full quantitative details are difficult to ex- 
tract. There are, however, some exact results for x\. the 
first maxima scales with ./V as (x±) ~ In AT and, further- 
more, the probability distribution P(y) for the shifted 
maximum y = X\ — (xi) is given by [38| 



P(y) = [2e y / 2 K 1 (2e y / 2 )]' 



(15) 



where K\ is the modified Bcsscl function. The agreement 
between P(y) and numerics in FigQJc) is excellent. 

As for the width of the distributions and the scaling 
of the gaps for a = 1, we only have indications from 
simulations. It appears that &k>\ ~ 0(1), as for k = 1, 
and that the gaps dk are also 0(1) for all examined k. 
Logarithmic crossovers prevent us from making any firm 
conclusions. 

For a > 1, the convergence to the limit distribution 
becomes power law as seen in extensive simulations for 
general a 24| and also in the exact solution for a = 2 
]iqL |23| . Thus the convergence is faster and the order 
statistics can be reliably investigated for a > 1 (but not 
too close to a = 1). 

FigQId) displays the a = 2 results for Pk(z), in which 
one can observe a new scaling regime. Namely, the his- 
tograms are barely distinguishable, i.e. the widths of the 
histograms a k are much larger than the gaps dk- It is 
known [l9| from the exact solution for P\{x) for a = 2 
that o\ ~ (xi) ~ N 1 / 2 . Thus our simulation results (see 
Fig© indicating (xi — x-i)ja\ ~ N^ 1 ^ 2 confirm that the 
first gap di = (xi — xq) is indeed 0(1). It then follows 
that there is no need for an iV-dependent rescaling of the 
gap. The same result is obtained for all dk — (xk ~ Xk+i) 
as seen directly in FigJU where the unsealed gaps for 
a = 2 are plotted as functions of N. 




N 



FIG. 4: (Color online) Unsealed gaps for a — 2. The simu- 
lation results are for N = 16,32, 16 384, and the lines are 
just guides for the eyes. 

Plots similar to those in Fig [J] implying that dk ~ 0(1) 
in the large N limit can be produced for any 1 < a < 
5. The same conclusion can also be drawn from Fig 12 
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Indeed, it is known [24| that o\ ~ N^ a x ^ 2 , and since 
the asymptotes of the simulation results for dk/&x scale 
like N-^- 1 !/ 2 , one concludes that d k ~ 0(1). 

Once dk ~ O(l) is established, we can plot the un- 
sealed dk as a function of k. As can be seen in Figj3l 
the dk display power law scaling and the exponents for 
a = 2, 3, 4 arc —1/2,0, and 1/2, respectively. These 
exponents, together with the —1 and +1 values (with 
possible logarithmic corrections) at a = 1 and 5, suggest 
that in the interval 1 < a < 5 

d k ~ k^-V' 2 . (16) 

The above scaling can be derived from the phenomeno- 
logical considerations embodied in cqs. (fT2"j) - (fl"4"l) which 
worked for a < 1. The starting point is the observa- 
tion that the sum of dk is of the order of x\, and since 
x\ ~ A^" -1 )/ 2 , we can write 

N/2 

Y^dk-i^-N^I 2 (17) 

fc=i 

Assuming that the dfe ~ A: 5 scaling is valid for large k of 
the order of N, we have 

N/2 

^Jfc'wJVC - 1 )/ 2 (18) 
fe=i 

and the above equality implies dk ~ fc( a_3 '/ 2 , as sug- 
gested by simulations. 

The large /c asymptote of the spectrum Sk can now 
be obtained by integrating the dk ~ fc( Q_3 )/ 2 expression, 
yielding 

£fc ~ k^' 2 . (19) 

As can be seen from Figj5l fits of the form Sk = 
a l t ia-i)/2_^f ) are m( Jeed excellent for a = 2, 3, 4. It should 
be noted, however, that the constant b is important for 
observing the exponent (a — l)/2. Naive straight-line 
fits on log-log plots of Ek versus k yield poor estimates 
for the exponents since the constant b generates strong 
corrections to scaling in the range < k < 10 studied. 

C. Large a regime (5 < a < oo) 

The large a regime is special in that the second deriva- 
tive of the signal becomes continuous for a > 5. Then the 
largest local maximum x\ = x(ti) and the values near it 
are expected to determine the order statistics, and this 
allows one to calculate the spectrum by expanding the 
signal around the maximum 

x(h + A) = xi + ^z"(fi)A 2 + 0(A 3 ) . (20) 

where A = r, 2r, (fc — l)r are the positions of the 
second, third, fcth maximum respectively. The justi- 
fication for the above expansion comes from comparing 
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k 



FIG. 5: (Color online) Spectra of the ordered values (a;*.) 
measured from the largest value of the signal Ek = (xi — Xk) 
for a = 2, 3, and 4. Simulation results for iV = 16 384 are 
plotted together with two-parameter fits (solid lines) of the 
form Ek = afc' cl_1 ' / ' 2 + b where a and b depend on a. 

the first and second terms on the right hand side of (|20p . 
The maximum value X\ scales as (xi 

) ~ iv^- 1 )/ 2 , and 

the characteristic value of x" can be estimated through 
[((a;") 2 )] 1 / 2 . A straightforward calculation yiels 

[((x") 2 )] 1/2 ~ (|x"(ti)|) ~ N^/ 2 ~ (xOA^- 2 (21) 

showing that the second term is indeed much smaller 
than the first one (note that A-r = T/A r -(9(l)in 
the limit we are considering). 

Using the expansion (|2"0|) . the unsealed gap dk can be 
written as 

{x{t x + (k— l)r) - x(ti + fcr)) w (\x" (ti)\)T 2 k . (22) 

Next, the iV-dependence contained in (\x"(ti)\)T 2 ~ 
N ( a -5)/2 jg sca i c d out, and we obtain the gap as a func- 
tion of k 

d k ~k (23) 

and the corresponding spectrum then follows as 

e k - k 2 . (24) 

As can be seen in Figj3l the a = 5 and 6 results for dk 
are in excellent agreement with the dk ~ k scaling for 
small values of k as well. 

It should be noted that the scaled first gap shown for 
a = 5 and 6 on Fig J2] can be also understood on the basis 
of the expansion (|20H22| . Indeed, noting that (xi) ~ ax, 
and using (f2Tj) with k = 1, we have 

(Xjtj-Xjh+T)) (\x"(h)\) 1 

ax ax 7V^' (25) 

in agreement with the a > 5 results on Figf2] The above 
result is exact in the a — > oo limit where the signal is 
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just a single mode x(t) ~ X\ sin (2nt/T), thus the rela- 
tionships (xi) ~ crj and (|x"(ti)|) ~ (x\)/N are trivially 
satisfied. 

The simulations for a = 5 and 6, together with the 
simplicity of the a — > oo limit, suggest that the results 
()23H25[) are valid for all a > 5. The possible logarithmic 
corrections at the borderline point a = 5 arc outside the 
limits of our present simulations. 

V. SUMMARY AND A CONNECTION TO THE 
DENSITY OF NEAR EXTREME STATES 

The results for the gap dk and the spectra Sk obtained 
in Sec lIVI are summarized in columns 1-3 of Table [U As 
one can see from columns 2 and 3, there are three regimes 
and, furthermore, the gaps and the spectrum change with 
the noise properties (a) only in the 1 < a < 5 regime. 

TABLE I: Gaps, dk and spectra Ek in the order statistics 
of l/f a signals as obtained in simulations and suggested 
by scaling arguments. The fourth column displays the one- 
dimensional potential V(x) in which the quantum mechanical 
motion of a particle generates energy spectra with the same 
scaling properties. 



a 


d k 




V(x) 


< a < 1 




Ink 


In \x\ 


1 < a < 5 


fc(«-3)/2 


fct"- 1 )/ 2 


i |2(ct-l)/(5-c0 


5 < a < oo 


k 


k 2 


\x\°° 



The above results for 1 < a < oo can be derived 
by considerations related to the density of near extreme 
states. Namely, if we assume that the spectrum has 
a scaling form (x\ — x k ) = £k ~ & then the density 
of states near the extreme is p(Sx) ~ (Sx) 1 ^ / 'Sx ~ 
(Sx) 1 ^^ 1 , thus knowledge of the small argument asymp- 
tote p(Sx) ~ (Sx) 1 allows one to deduce the exponent 

= 1/(7+1)- 

The density of near extreme states p(Sx) was first in- 
vestigated for i.i.d. variables by Sabhapandit and Ma- 
iumdar (2o| . Later, it was shown by Burkhardt et al. 
j26| that p(5x) can be obtained for periodic signals as 
the distribution $i(x) of the maximum with respect to 
the initial value x = max t a;(£) — x(0). The distribution 
&i(x) for periodic 1//" signals has been investigated by 
simulations as well as through exact solutions for the par- 
ticular cases of a = 2, 4, and oo (26|. The small argument 
behavior of was indeed found to have a power law form 
$z(5x) - (fa) 7 with [see eq.(65) in [H] ] 

— — - , a < 5 , 
1 (a)=l a i (26) 

1-2' 



The scaling exponents of the spectra following from the 
above expression [f3 = 1/(7 + f)] are equal to those in 
Table U for 1 < a < 00, thus reinforcing our earlier con- 
clusions. Since the exponents in (|26p are exact for a = 2, 
4, and 00, we can conclude that the corresponding (3(a) 
exponents depend on the validity of the e k ~ scaling. 

VI. COMPARING WITH QUANTUM SPECTRA 

We shall now compare the order statistics spectra to 
the energy spectra of quantum mechanical systems in 
the quasi-classical limit. The reason for this comparison, 
apart from its entertaining aspects, is that the discrete 
quantum mechanical spectra may also be considered as 
an order statistics spectra. 

Let us consider a particle of mass m which moves in a 
potential 

V(x)=g\x\ 6 (27) 

where g > is the coupling constant. The simplest way 
to calculate the quasi-classical limit of the spectra is to 
use dimensional analysis combined with the observation 
that, in the large quantum-number limit (k — > 00), the 
quantization condition (Jpdq = kh) forces h and k to 
appear in the combination hk. This means that the k- 
dependence of the spectra is determined by its h depen- 
dence. Since the energy is uniquely determined by the 
dimensions of m, q, and h, one obtains 

E k ~ (hk)-&* . (28) 

Comparing the above result with Column 3 in Table [H 
one can see that there is a mapping between the quantum 
mechanical and the order statistics exponents. The large- 
a regime (5 < a < 00) corresponds to the infinite square- 
well potential (9 = 00), while 9 changes monotonically 
from to 00 in the nontrivial regime (1 < a < 5), and 
the correspondence is given by 

6 = ^1. (29) 

5 — a 

The above considerations cannot be applied for the re- 
maining < a < 1 range but the order statistics spec- 
trum £fc ~ In k suggests that the corresponding potential 
is of the form V(x) ~ In 

We would like to emphasize that it is not obvious that 
the mapping between the exponents have any physical 
content. Nevertheless, it is intriguing to ask whether 
there is a quasi- classical extreme value question whose 
answer is the quantum mechanical spectra. 

VII. FINAL REMARKS 

Studying the order statistics in 1// Q signals, we found 
three scaling regimes as summarized in Table |TJ In the 
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case of weakly correlated stationary signals (0 < a < 1), 
the i.i.d. result (dk ~ k~ ) applies, while for strongly 
correlated signals (1 < a < oo) the scaling depends on 
whether the signal is twice differcntiable (a > 5) or not 
(1 < a < 5). In the former case, order statistics (dk ~ &) 
follows from expanding the signal around the maximum. 
In the latter case, the observed scaling (dk ~ fc^ Q_3 ^ 2 ) 
is derived using a phenomcnological argument (exact for 
i.i.d. variables), namely that the sum over dk scales in 
the same way as (x\). Meanwhile, the spectrum Sk is 
essentially obtained by integrating d^. The same scaling 
picture was also obtained by relating order statistics to 
the density of near extreme states studied in previous 
works. 

It is clear that investigating the order statistics of cor- 
related systems will help in characterizing and under- 
standing extreme events in more detail. It is also clear, 
however, that the results for l/f a signals have a re- 
stricted range of applicability, and much more is needed 



to advance our understanding of the effects of various 
correlations. Nevertheless, there are phenomena where 
1// Q type fluctuations do emerge, the extremes are im- 
portant and, consequently, our results may be utilized. A 
possible application is to climate records (temperature, 
precipitation, etc.) which often display e.g. clustering of 
extreme events, and the power spectrum of the fluctua- 
tions appears to be of 1// Q type [1, S3]. 
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